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ABSTRACT 

We consider the issue of characterizing the coherent large-scale patterns from CMB 
temperature maps in globally anisotropic cosmologies. The methods we investigate are 
reasonably general; the particular models we test them on are the homogeneous but 
anisotropic relativistic cosmologies described by the Bianchi classification. Although 
the temperature variations produced in these models are not stochastic, they give 
rise to a "non-Gaussian" distribution of temperature fiuctuations over the sky that 
is a partial diagnostic of the model. We explore two methods for quantifying non- 
Gaussian and/or non-stationary fiuctuation fields in order to see how they respond to 
the Bianchi models. We first investigate the behavior of phase correlations between the 
spherical harmonic modes of the maps. Then we examine the behavior of the multipole 
vectors of the temperature distribution which, though defined in harmonic space, can 
indicate the presence of a preferred direction in real space, i.e. on the 2-sphere. These 
methods give extremely clear signals of the presence of anisotropy when applied to 
the models we discuss, suggesting that they have some promise as diagnostics of the 
presence of global asymmetry in the Universe. 
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1 INTRODUCTION 



Observations of the Cosmic Microwave Background (CMB) provide some of the most compelling support for the cur- 
rently favored ACDM, or concordance^ cosmological model. The concordance framework predicts that the CMB should 
posses temperature fluctuations which a r e both statistically isotropic (i.e. st ationary over the celestial sphere) and Gaus- 
sian (jGuth Pilll982l : Istarobin skiil 1982 : Bardeen, Steinhardt Turner 1983l V Measurements by the Wilkinson Microwave 
Anisotropy Probe (WMAP) (Ben nett etaD 1 20031 : Einshaw et al.ll2009l ) have under gone extensive statistical analysis, much of 
which has confirm e d the concordance model but with some indications of departures that may be significa nt; see for example 
Yadav Wandeltl (I 20081) . More specifically, there is some evidence fo r hemispherical power asymmetry (jEriksen et al.l [20041: 



Parlj|2004l:[Eriksen et al.ll2007l : Eoftuft et al.ll2009l : iHansen et al.ll2009l ) and also a Cold Spot has been identified (|Vielva et al 



20oi ICruz et al.ll2005l ). In other words there is some evidence of an anisotropic universe, i.e. one in which the background 



cosmology may not be described by the standard Friedman-Robertson- Walker (FRW) metric. Of course the background cos- 
mology for a non-isotropic universe may still be described by the FRW metric, but this would require a non-standard topology 
which we do not consider in this analysis. 

The Bianchi classification provides a complete characterization of all the known homogeneous but anisotropic exact 
solutions to General Relativity. The classification was first proposed by Bianchi and later applied to General Relativity 

( Ellis M acCallum 1969). Initial stud ies used the lack o f large-scale asym metry in the CMB temperature to put strong con- 

straints on the possible Bianchi models ( Barrow. Juszkiewicz Sonoda 19851 : Bunn. Ferreira Silk 1996 : Kogut. Hinshaw Bandav 



19971 ). However, simulations of the CMB from Bianchi universes not only show a pr eferred direction, but models with neg ative 
spatial curvature (such as the types V and VII^) can produce localized features ( Barrow. Juszkiewicz Sonodal 1985 ). So 
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more recently attention has shifted to reproducing a Cold Sp ot such as that claimed to exist in the WMAP data. Initially, Type 
Vllh, was the favored mod el to best reproduce t he anomalv (|jaffe et al . 2005, 2006a fl), a nd thi s has subsequently been inves- 



Pontzen Challinorl I2OO7I : IPontzenI lioool : 



tigated quite thoroughly (|McEwen et al. I I2QQ6I : [Bridges et al.ll2008l : 
2010) , although more recent work has also looked at the Bianchi Type V which also produces localized features 
2OO9II. 
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The most interesting range of anisotropic structures is produced in Bianchi Types VII^,, VIIo and V. These different 
Bianchi types have the effect of focusing and/or twisting the initial quadrupole over time (see Figure [Tj). In this paper we 
study the behavior of these Bianchi models so as to identify characteristics of the radiation fields they produce and develop 
methods that can be used to identify more general forms of anisotropy. Understanding the characteristics identified in these 
particular cases will hopefully help us find better and more systematic ways of constraining the level of global symmetry 
present in the real Universe. Note we consider just characteristics observable in the CMB temperature; we shall return to a 
study of the polarization radiation component in later work. 

We consider two statistical measures of anisotropy in some detail in this paper. Neither of these is entirely new and both 
have previously been applied to observed CMB maps. However, the general philosophy behind previous applications of these 
methods has been simply to look for departures from the (composite) null hypothesis of statistica l isotropy and Gaussianity ( or 
more recently they have been developed to look at universes with multiply-connected topologies (|Bielewicz k Riazuelolbood ll. 
In other words, they have been used to construct hypothesis tests with the concordance cosmology but their performance has 
not hitherto been evaluated on models with built-in anisotropy. 

For example, if the concordance model is correct, the phases of the spherical harmonic coefficients of t he CMB should be 
independently random and uniformly distributed. Recent studies have suggested some deviation frora this (jColes et al.ll2004l : 
Stannard k Coleslbooj : IPineen. Rocha k Colesll2005l : Ichiang. Naselskv Colejl2007l : IChiang et al.ll2007l ) but it is not clear 



whether they indicate global anisotropy or departures from Gaussianity, let alone whether these are of cosmic or instrumental 
origin. Here we examine the use of phase correlations in quantifying the temperature patterns generated in models with known 
levels of global inhomogeneity. 

There have since bee n attempts to under- 



Multipole vectors were first introduced over a century ago ()Maxwelllll891l ) ^ 

stand the multipole vectors in order to explain the CMB anomalies reported a t large angular scales (|Katz k Weeksl 20041 : 
Schwarz et all I2OO4I : ICopi. Huterer StarkmanI I2OO4I : iLand Magueiiol boOSal lblHIdlH: ICopi et al.1 12OO6I . I2OO7I ) since is not 



clear how to quantify and verify such properties from the CMB anomali es in spherical harm onics. They have been used in 
a num ber of studies to show anomalie s, such as alignments of multiples ( Abramo et al. 20061 ) in a similar plane to the axis 
of evil (jLand Magueiiol l2005dl . I2OO7I ). Our aim here is to examine the behavior of the multipole vectors in cases where the 
form of anisotropy is known priori in order to assess their potential to act as more general descriptors. 

Two points are worth making before we continue. First, any realistic cosmology (whether of FRW or Bianchi type) 
will possess random fluctuations on top of a smooth background. If these fluctuations are stationary Gaussian then they 
will add correlated "noise" to any signal arising from the background model and will thus hamper the performance of any 
statistical analysis method, especially at smaller angular scales. This Gaussian "noise" (which is equivalent to stationary 
Gaussian fluctuations, and not to be confused with instrumental noise) is completely characterized by second-order statistical 
quantities (i.e. the power spectrum in harmonic space or the autocorrelation function in pixel space). The statistical descriptors 
we explore are independent of the power-spectrum, so adding Gaussian noise will not produce any systematic response in them. 

We also restrict ourselves to looking at just the large-scale features because the patterns in the temperature maps 
resulting from the Bianchi models is over large scales. Therefore, by looking at large scales only, there is more chance of 
detecting the anisotropy. However, it goes without saying we are not claiming that these Bianchi models are in themselves 
complete alternatives to the concordance cosmology. Rather we think of them as representing possible perturbation modes of 
the FRW background. 

The layout of this article is as follows. In Section [2] we look at pixel distributions of the CMB maps to show how the 
statistical anisotropy present in these models produces a form of non- Gaussianity in the pixel distribution over the celestial 
sphere. We then introduce phase correlations in Section [3] to provide characterization of the anisotropy displayed by the 
models. In Section |4] we look at the behavior of the multipole vectors as characteristics of the anisotropy of the same maps. 
Finally, Section [5] summarizes the conclusions. 



2 STATISTICAL DESCRIPTORS OF NON-GAUSSIANITY 



Sung ColesI ((2009) discussed how localized features in the CMB temperature pattern, perhaps similar to the Cold Spot 



observed in the WMAP data dvielva et al.ll2004l : ICruz et al ] l2005h . can be generated in models with negative spatial curvature, 
i.e. Bianchi types V and VII/^. In the standard cosmological framework the temperature fluctuations are described by a 
Gaussian random process over the sky, so a feature like the cold spot corresponds to an extreme event in the tail of the 
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Figure 1. Simulated maps of the the CMB temperature, at redshift z = 0, using Bianchi type cosmologies. From left to right the Bianchi 
types are: V, VIIq and VII/^. The colour scale is marked in milliKelvin. All the maps started as a quadruple at z = 500. The Bianchi V 
map shows a focused feature, the Bianchi VIIq map has a twisted feature and the Bianchi VII/j, map has both focusing and twisting in 
the resulting temperature pattern. 

distribution of fluctuations. In a Bianchi model, however, it is not stochastic at all but produced coherently as a result of the 
geometry of the space-time. 

Clearly we need a systematic way to characterize the relationship between rare events like this and their origin through 
either non-Gaussianity or global anisotropy. Analysis of the temperature patterns using standard descriptors in non-standard 
scenarios will produce signals different from what one would see in the presence of stationary Gaussian noise. To illustrate 
this issue we study the pixel distribution function of temperature values. This is an obvious way to test for non-Gaussianity 
in a random field of temperature values, but a coherent fluctuation field also possesses a one-point distribution that yields 
some diagnostic information. In this sense, all the Bianchi models are inherently non-Gaussian but their non-Gaussianity is 
simply a manifestation of the presence of anisotropy. 

We calculated the pixel distribution function, which is simply a frequency count of the pixel (or temperature) values, for 
each of the Bianchi maps and used it to plot the histogram seen in Figure [2] A perfectly homogeneous and isotropic map, such 
as that predicted by the concordance model, would have a constant value over the whole map - remember we are considering 
maps without fluctuations - and therefore a histogram of the pixel distribution histogram for this map would give a delta 
function at the mean. Our results show some deviation from this prediction. In Figure [2] we see the plots for the Bianchi V 
and VII^ types have strongly peaked features at the mean, but still with some non zero variance; the type VIIo model has a 
nearly uniform distribution across the whole temperature range. 

Although not demonstrated in this diagram, another point to note is that at early times the histograms of the pixel 
distribution functions of the three different Bianchi types are almost identical; the different values of the temperature pixels 
are roughly even over the range. As redshift^ decreases, the temperature patterns for Bianchi V and VII^ gradually start 
to focus fsee lSung Colesll2009l ). and their histograms of the pixel distribution functions become successively more peaked. 
For the Bianchi type VIIo the temperature pattern just twists, reorganizing the pattern on the sky while the histogram stays 
roughly the same. These observations help to explain the features in the histograms. As the temperature patterns become 
more focused, more of the rest of the map becomes uniform and so the histogram is tending towards a delta function. 

In summary, what we would hope to discover from the pixel distribution histogram is that it gives us some clues about the 
homogeneity and isotropy of the maps i.e. the degree of concentration around the mean might tells us about the homogeneity 
of the parameters or the asymmetry of the distribution might give information about the anisotropy. But as it stands, the 
information from the pixel distribution function is not that clear. All we can say is that the histograms differ from a delta 
function, so the maps are not perfectly homogeneous and isotropic, and that the histograms are clearly non- Gaussian in shape. 
However the shape of this one-point pixel distribution does not furnish us with a complete description of the pattern because 
it does not take into angular correlations between the pixels. 



3 PHASE CORRELATIONS OF BIANCHI CMB MAPS 

We now move on from pixel distributions to consider the spherical harmonics of the temperature maps, or more specifically 
the phases of the spherical harmonic coefficients. 

3.1 Spherical Harmonics 

The temperature of the CMB, Te),^, is defined on a sphere where ^ [0? ^] and G [0, 27r] are the polar and azimuthal angles. 
Therefore one way of describing the temperature anisotropics, ATe),^, is to extract the corresponding spherical harmonic 

^ Note, redshift is defined here as proportional to the inverse geometric mean of three scale factors. 
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Figure 2. The pixel distribution function of the normahzed temperature for the Bianchi CMB maps. The y axis indicates the normahzed 
pixel number which is the whole range of temperature covered by the plot, but normalized so they can be plotted on the same axes. The 
colour cyan represents type V, magenta for VII/^, and green for VIIq. The plot shows that the Bianchi V and VII/^ types have strongly 
peaked features at the mean, whereas the type VIIq has a nearly uniform distribution across the temperature range. 



coefficients (a^m): 
Te,ct) — T 



AT, 



(1) 



£=0 m= 



where \airn \ and are the amplitudes and phases of the spherical harmonic coefficients, and Y^rn are the spherical harmonics 
which are defined here as: 



21 + 1 {l-m)\ 
47r {l + m)\ 



(2) 



where Pira is the associated Legendre Polynomial. Note that this definition of spherical harmonics includes a phase factor of 
( — 1)"^, also known as the Condon-Shortley phase. 

In the standard cosmological model, the temperature fluctuation field is produced by stochastic fluctuations which are 
Gaussian and statistically stationary over the celestial sphere. In t his case the phase s ^ira of each spherical harmonic mode 
a^rn are independent and uniformly random on the interval [0, 2ii] (|Coles et al. 12004 V If instead the temperature pattern on 
the sky is produced by a Bianchi geometry then the airn are no longer stochastically generated but can be directly calculated 
from pa rameters of the model. Analytical forms for the temperature pattern can be used to obtain the spherical harmonic 
phases ([McEwen et al. 20061 : Bridges et al. 20081 ). but it is clumsy to transform these between different coordinate systems 
( Coles et al.l 2004ll . In the followin g we therefore obtain distributions of from Bianchi maps generated using the method 
described bv lSung Coled (|2010l ). 



3.2 Visualising phase correlations 

To visualize the information held in the phases, ^^m, of the spherical harmonic coefficie nts, a^ra, we plotted th em over all 



and m. Rather than using a 3D plot, colour has been used to represent the following Coles Chiang! ( 2000l ) . The colours 



equate to the angle on a colour wheel: red (^^m = 0), green (^^m = ^/2), cyan (^^m = tt), and purple {^im — 37r/2). To 
understand these plots, first consider what we would expect to see in the case of an isotropic and homogeneous universe as 
predicted by the concordance model. This would be a uniform map (as we are not at this point considering fluctuations) but 
in spherical harmonics this only has power in one mode = m = 0), so there is no phase for the other modes. Better to 
consider a map with Gaussian fluctuations as later in the section we will move on to add noise to the Bianchi maps. Figure 
[3] shows the phases ($£m) for a homogeneous and isotropic map with Gaussian fluctuations. The phases are random over the 
space i.e. there are no visible patterns in the distribution of colours in the plot. Note that for all the maps, ^^rn = or tt for 
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Figure 3. Example of the spherical harmonic phases (^£rn) we would expect for the concordance model where m ^ [0, 20]. The 
distribution of ^£rn is random, apart from for m = 0, where the phases can only be or tt by definition. The colours represent the 
different values of the ^£rn (red (^£rn = 0)? green (^£rn = 7r/2), cyan (^£rn = ^r), and purple (^£rn = 37r/2)). 




Figure 4. Phases of the spherical harmonic coefficients (^im) for Bianchi types V (left), VIIq (middle) and Yllh (right) where m G 
[0, 20] and z = 0. Note that £ is plotted against the x axis, increasing from left to right, and m is plotted against the y axis, increasing 
from bottom to top. The distributions are not random (as in Figure [Sj but exhibit some distinctive features. All the ^irn for the VIIq 
and V types are orthogonal, and there are sequences of colours in the type VII^ (see m = 2). 



Figure 5. A^^^n for Bianchi types V (left), VIIq (middle) and Vllh (right) where m G [0, 20] and z = 0. Note that i is plotted against 
the X axis, increasing from left to right, and m is plotted against the y axis, increasing from bottom to top. Like the phases (Figure 3]) 
the distributions are not random but exhibit some distinctive features. All the A^£^ for the V type are either or tt. The A^£^ for 
the VIIq type are again orthogonal but in a more correlated way. Similarly, the sequences of colours in the type VII/^ are now even more 
prominent (see m = 2). 



m = because the aim coefficients are defined so that aim — ai^-m- Other than this, the distribution of is random. Also, 

note that for all \m\ > — 0. 

We extracted the from each of the Bianchi maps using Healpix (jGorski et al.ll2005l ) and plotted them in the same 
way as we have described; the results are shown in Figured] The plots show that the are not random but have patterns, 
i.e. the harmonic modes manifest some form of phase correlation. For all the Bianchi types, — for all odd m. For the 
VIIo and V types, all the ^irn are orthogonal i.e. they are either 0, 7r/2, tt, or 37r/2. Both the VIIo and VII^ types show 
sequences of increasing/decreasing phases, which are particularly prominent for m = 2. 
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Table 1. Results from the Kolmogorov-Smirnov test comparing the distribution of phase differences in the Bianchi CMB maps with a 
random distribution of phases as predicted by the concordance model. D is the Kolmogorov-Smirnov statistic found by comparing the 
phase differences (A^). P{D) is the Monte Carlo estimate of the probability of getting the value of D, or less, found for the Bianchi 
models, from a random selection of phase differences. These are computed by forming an empirical distribution of D from sets of random 
simulations and counting what fraction of the ensemble gives the results obtained for the Bianchi maps. For example, in the case of the 
P{D) for the VII/, map (z = 500) we find that, out of 10000 simulations, 9450 have a value of D less than 0.11. Given the probable 
sampling accuracy of around one percent, we have rounded the results. 

While some patterns are apparent in these plots, an even better way to visualize the phase correlations is to look at the 
phase differences which are defined here as: 

A^im ^im - ^i,m-l (3) 

The phase differences are shown in Figure [5] and the correlations are much more apparent compared to the plots of ^im- All 
the A^im for the V type are lined up, i.e. either or tt. The A^im for the VIIo type are again orthogonal, but whereas in 
the phases the distribution of 0, 7r/2, tt, and 37r/2 seemed some what random, in the phase differences we see similar values 
'clump' together. Similarly, the sequences of colours in the type VII/, (see m = 2 for example) are now even more prominent. 

So, strong correlations are observed in the phases and phases differences of the simulated Bianchi CMB maps. But we 
have only looked at large angular scales where there are only a small number of independent data points. Even without noise, 
it is important to ask the question whether these correlations are likely to be statistically significant. One way to quantify 
this is to use a Kolmogorov-Smirnov test. This is a non-parametric statistical test which measures the maximum distance of 
a given distribution from a reference probability distribution. In this case we want to show deviation from a random set of 
A$£rn, i.e. a uniform distribution, which is predicted by the concordance model. 

To calculate the Kolmogorov-Smirnov test statistic a set of phase differences A^^^n are separated into bins of equal size 
between and 27r. The number of A^im which fall into each bin are counted and a cumulative distribution derived. If the 
distribution is uniform, as in the case of the reference probability distribution, then the number of A^im in each of the bins 
should increase roughly linearly across the bins. The difference between both the sample and uniform cumulative distributions 
is found for each bin and the biggest difference is the Kolmogorov-Smirnov statistic D. 

To deduce the significance of Z), a set of ten thousand tests have been run to generate sets of random angles of equal 
size to the sample sets. D was found for each of these sets and this data was used to find the significance of D for the sample 
distributions from the Bianchi maps. 

The Kolmogorov-Smirnov statistic D, and the derived probability of that statistic P{D), for all the Bianchi maps are 
detailed in Table [1] This table shows that there is indeed a significant deviation from a uniform distribution for the phase 
differences for all Bianchi types. Of the 10000 random sets of data, none showed a value for D as high as seen for the Bianchi 
cases. 

The Bianchi VII/, type was also considered at different redshifts to see how the correlations changed with time. Table [1] 
shows that in general value of D gets more significant over time i.e. the correlations in the phase differences of the Bianchi 
maps become stronger over time. 

3.3 Rotating maps and adding noise 

In Section [3. 2 1 we applied the Kolmogorov-Smirnov test to a "clean" map that is perfectly aligned with the vertical axis. This 
section addresses how noise and rotation affect the identification of correlations in the phases of the spherical harmonics of 
CMB maps from Bianchi models. 

First we consider rotation. Phases of spherical harmonic coefficients are not rot at ion- invariant. Rotating the coordinate 
system used to represent a CMB map in (f) (which is equivalent to rotation around the z axis) would increment each of the 
spherical harmonic phases by 0, so the phase differences would remain the same. Therefore rotation in (f) would have no effect 
on the value of the Kolmogorov-Smirnov statistic D. Rotation in 6 is more complicated to express so we used an empirical 
approach to quantity the effect on D. The Bianchi CMB maps were rotated by a small angle, 6 — 7r/8, and then the spherical 
harmonic coefficients were derived and used to calculate D. The results in Table [2] show that the values of D for each of 
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Table 2. Results from the Kolmogorov-Smirnov test comparing the distribution of phase differences in the Bianchi CMB maps, rotated 
by = 7r/8, with a random distribution of phase differences as predicted by the concordance model. D is the Kolmogorov-Smirnov 
statistic found when considering the phase differences. P(D) is the Monte Carlo estimate of the probability of getting the value of -D, or 
less, found for the Bianchi models, from a random selection of phase differences. These are computed by forming an empirical distribution 
of D from sets of random simulations and counting what fraction of the ensemble gives the results obtained for the Bianchi maps. For 
example, in the case of the P(D) for the Yllh map we find that, out of 10000 simulations, over 9999 have a value of D less than 0.46. 




Figure 6. A^i^n for m e [0, 20] (Bianchi type V map at z = with white (left), WMAP (middle), and ACDM fluctuations (right) 
noise maps , = 7r/8). Note that i is plotted against the x axis, increasing from left to right, and m is plotted against the y axis, 
increasing from bottom to top. Correlations can still be observed as lines of similar colours. 



the maps are even higher than in maps that hadn't been rotated, indicating the presence of even stronger correlations. This 
suggests that, at least for small rotations off the axis, the correlations are just as significant, if not more so. 

As an aside, the colour plots of the phase differences for Bianchi maps rotated by a number of different in the range 
to 27T were generated. These plots have been condensed together into movies^ which show that the correlations in the VIIo 
and V maps are visible across all and for the VII^ map are visible within about 7r/3 of the preferred axis. . 

Now to investigate the effect of noise, we considered three different types of noise. Firstly we tried the simplest form by 
just adding white noise to the Bianchi map. A map of random Gaussian noise (white noise) was generated. Using Healpix 
the spherical mode resolution was reduced to ^ ^ 20. Then the "noise" map was modified to have zero mean and variance 
half that of the Bianchi map. The second "noise" map was derived from a product available on the WMAP Lambda^ website 
which provides the effective number of observations per pixel. A map of random Gaussian noise was again generated. The 
variance was modified per pixel so that it was inversely proportional to the square of the number of observations in that pixel. 
Using Healpix the spherical mode resolution was reduced to ^ ^ 20. Then the noise map was modified to have zero mean 
and variance half that of th e Bianchi map. The final "noise" map used a simulation of ACDM fluctuations of the CMB (as 
performed bv lSriksen et al. Again the noise map was modified to reduce the spherical mode resolution to ^ ^ 20 and 

have variance half that of the Bianchi map. 

Each of these "noise" maps was added to each of the rotated Bianchi maps. We see from the example in Figure [6] that the 
spherical harmonic coefficients derived still have visible correlations in the phases for the Bianchi V map. The results of the 



^ The movies can be found at |http: / /www. astro. cardiff.ac.uk/research/theoreticalcosmology/?page=research| 
^ http://lambda.gsfc.nasa.gov/ 
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Table 3. Results from the Kolmogorov-Smirnov test comparing the distribution of phase differences in the Bianchi CMB maps rotated by 
6 = 7r/8 with white, WMAP, and ACDM noise maps (z = 0) . D is the Kolmogorov-Smirnov statistic found when considering the phase 
differences. P{D) is the Monte Carlo estimate of the probability of getting the value of D, or less, found for the Bianchi models, from 
a random selection of phase differences. These are computed by forming an empirical distribution of D from sets of random simulations 
and counting what fraction of the ensemble gives the results obtained for the Bianchi maps. For example, in the case of the P(D) for 
the VII^ map with white noise we find that, out of 10000 simulations, 9920 have a value of D less than 0.15. 
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Table 4. Results from the Kolmogorov-Smirnov test comparing the distribution of phase differences in the WMAP ILC map, rotated 
to align with either the galactic axis or axis of evil, with a random distribution of phase differences as predicted by the concordance 
model. D is the Kolmogorov-Smirnov statistic found when considering the phase differences. P{D) is the Monte Carlo estimate of the 
probability of getting the value of or less, found for the Bianchi models, from a random selection of phase differences. These are 
computed by forming an empirical distribution of D from sets of random simulations and counting what fraction of the ensemble gives 
the results obtained for the Bianchi maps. For example, in the case of the P(-D) for the ILC map in the galactic plane we find that, out 
of 10000 simulations, 6620 have a value of D less than 0.06. 



Kolmogorov-Smirnov test (see Table [3]) show that the correlations are still detectable and significant for the Bianchi V and 
Yllh maps but not so well for the VIIo maps. So the method is better for detecting focused features that twisted features. 

We see that the effect of adding fluctuations here is not dissimilar to adding just Gaussian noise. The concordance 
model predicts fluctuations which are stationary and Gaussian, as discussed in the Introduction (Section [1]). Although these 
fluctuations are correlated on the sky, they have random phases so are incoherent with respect to what our statistic measures. 

The "noise", or fluctuation, maps are added to the Bianchi maps so that the ratio of the variances is of order unity. 
However, any ratio is possible; this specific choice is just for illustrative purposes to demonstrate the proposed methods. 
Nevertheless, if a random-phase (Gaussian) signal is superimposed on the Bianchi template, the phase coherence of the 
resulting map will still be degraded. If the Gaussian component is too large, the overall map will be indistinguishable from 
one with purely random phases. 

In the examples we have shown, the Gaussian "noise" or fluctuation maps are added to the Bianchi maps in such a way 
that the ratio of the overall variance is of order unity. Our method still functions well with this level of "contamination" , but 
if the noise variance is much higher than that of the Bianchi maps the method begins to struggle. 

So summarizing, the phases of the spherical harmonic coefficients are a very effective way of identifying focusing features 
in CMB maps, as long as the noise is not excessive, and can be used to give quantifiable significances. However, like the pixel 
distributions, the variation from the expectation of the concordance model only gives us an indication of non-Gaussianity. It 
is not clear in what form the non-Gaussianity occurs, such as an anisotropy or inhomogeneity. Hence the next section looks 
at multipole vectors which are built from spherical harmonic coefficients but can be used to give results in pixel (as opposed 
to harmonic) space, which is more meaningful from the point of view of diagnosing the presence of a preferred direction. 



3.4 Application to WMAP 5 Year Data 

For pedagogical interest, the methods described in Section [321 are applied here to the WMAP 5 year Internal Linear Combi- 
nation (ILC) map. The results of the Kolmogorov-Smirnov test on the ILC map in the galactic coordinate system show very 
low significance correlations in A$ (see Table [4j). 

However Section [X3l showed that to see the phase correlations, the Bianchi maps needed to be rotated close to the preferred 
axis. There have been studies that have found a preferred direction in the WMAP data, highlighted by the alignment of at 
least the quadrupole {I = 2) and octopole (^ = 3). This preferred axis is known as the Axis of Evil. Therefore the methods 
from Section [3.21 are applied to the ILC map rotated so that the vertical axis aligns with the Axis of Evil. These A$ plotted 
in Figure [71 do not show any visual correlations. For comparison the figure also includes a plot of the same ILC data but 
with the phases replaced with random angles (i.e. so as to not affect the magnitude of the amplitudes of the ai^m)- Whilst 
the Kolmogorov-Smirnov test (Table [4]) finds higher significance results than when the map was in Galactic coordinates, the 
results are still at a low significance. 

The results show there is no significant detection, so if we do live in an anisotropic universe then it must be obscured 
with considerable "noise" (fluctuations). However the fact that the significance of the results do increase when the map is 
aligned with the Axis of Evil is intriguing. 



4 MULTIPOLE VECTORS FROM BIANCHI UNIVERSE 

As shown in the previous section, the properties of spherical harmonic coefficients provide us with a generally effective way 
of identifying anisotropy through correlations in CMB maps. However the geometric interpretation of the mode correlations 
seen in harmonic space is by no means easy to interpret geometrically. In an effort to use the spherical harmonics to provide 
more meaningful explanation of non-Gaussianities found, we now consider an alternative approach, based on multipole vectors 
These can be constructed from spherical harmonics, using the a^rn coefficients derived from CMB maps, but they give results 
in real (i.e. pixel) space. For a sum mary of the main term inology for the multipole vectors, using the polynomial interpretation 
approach which was introduced bv lKatz Weeksl (|2004l ). please see Appendix [Al 
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Figure 7. A$ for ^, m G [0, 20] for the ILC map with the Axis of Evil ahgned with the preferred axis (left) and for the same map but 
with the phases replaced with random phases (right). Note that £ is plotted against the x axis, increasing from left to right, and m is 
plotted against the y axis, increasing from bottom to top. No correlations are visible in either plot. 




Figure 8. The multipole vectors from Bianchi V map: dipole (left), quadrupole (middle) and octopole (right). 
4.1 Results for Bianchi maps 

Figure [8] shows the multipole vectors from the Bianchi V map which serves as a good example to show how strongly the 
multipoles are correlated. The dipole (left) lies exactly at the top of the sphere, which is at the centre of the image. The two 
quadrupole vectors (middle) are located on the same spots on which two of octopole vectors (right) are placed. The remaining 
octopole vector is in the centre, i.e. the same place as dipole. Now we plot the dipole, quadrupole and octopole on the same 
2-sphere for all the Bianchi maps, at different redshifts, to see how exactly they overlap (see Figure [9]). The background colours 
indicate if any of the multipoles over-lap; yellow for no overlapped multipoles, green for overlapped dipole and octopole, light 
yellow for overlapped quadrupole and octopole, and light blue if all the multipole vectors are overlapped. The z-axis is into 
the page and the x-y plane is the large marked circle. 




Figure 9. The multipole vectors from the Bianchi V (left), VII^ (middle) and VIIq (right) maps from early stage (bottom panel) to late 
time (top panel). The vectors are represented by dots: dark blue for the dipole, light red for the quadrupole and brown for the octopole. 
Background colours also indicate if any of the multipoles over-lap; yellow for no overlapped multipoles (right column and bottom of 
middle column), green for overlapped dipole and octopole, light yellow (bottom of left column) for quadrupole and octopole and light 
blue (top of left column) if all the multipole vectors are overlapped. 
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Figure 10. The multipole vectors from the WMAP 5 year ILC map: dipole (left), quadrupole (middle) and octopole (right). The top 
row shows results where the z-axis is into the page and the x-y plane is the large marked circle. The bottom row shows results where the 
z-axis is the vertical line and the x-y plane is the horizontal line across the centre of the map. 



First of all, in all types of the Bianchi models, we see the quadrupole and octopole vectors lie on the same plane, except 
for one of the octopole vectors which is always located in the centre of the image. For the Bianchi V and VII^, types, the 
dipole vectors lie very near the centre in the early stage but not exactly on it. However, as time goes on, the dipole vector 
is overlapped by one of the octopole vectors in the centre. The dipole vector of Bianchi VIIo type is different from Bianchi 
VII^. In the Bianchi VIIo type, there is no particular correlation between the dipole and other multipoles since the dipole is 
not coupled with the other multipoles (quadrupole and octopole). For the Bianchi V (left), at the beginning of the stage, the 
dipole vector is 'almost' on the z-axis and one of octopole vectors is exactly on the z-axis. However the two later-time Bianchi 
V cases on the top left of Figure [9] show the dipole and octopole are on the same spot, in the middle of the image, which 
is exactly on the z-axis. Meanwhile, one of components of quadrupole and octopole vectors are on the same spot on the x-y 
plan. This means that the later Bianchi V models have more extreme correlation between the multipoles. 



4.2 Results for WMAP Data 

Again for pedagogical interest, the multipole vectors are applied to the WMAP 5 year Internal Linear Combination (ILC) 
map in Figure [TOl At first glance, the results in the top row look like they might be clustered in a similar way to those in the 
Bianchi models (see Figure [9]). However when we look at the results in the second row, where the map is orientated in the 
usual galactic coordinates, we see the multipoles line up along the x-y plane. This suggests an alternative explanation that 
the clustering near the x-y plane could be a feature of the residual galactic contamination which is known to remain in the 
full sky ILC maps. 



5 DISCUSSION AND CONCLUSIONS 

The aim of this article was to explore some simple ways of characterizing the large-scale temperature patterns in CMB maps 
generated in anisotropic Bianchi type V, VII^ and VIIo universes. The ultimate purpose of investigating this behavior is to 
find ways of quantifying the global properties of the pattern produced in order to isolate the effect of anisotropy from that 
of non-Gaussianity. We repeat that when we talk about non-Gaussianity here is not related to a stochastic field; there are no 
fluctuations in the Bianchi maps. 

We first discussed perhaps the simplest and perhaps the most obvious possible descriptive statistics, the histogram of the 
pixel values, primarily with the aim of demonstrating how non-Gaussianity of a sort can arise from asymmetry. We evaluated 
the pixel distribution functions for each of the maps and compared them to results expected in a universe consistent with the 
concordance model. The type VIIo maps show the strongest deviation from the null hypothesis; but types V and VII^ behaved 
in a similar fashion to each other, and closer to that of the null hypothesis. The reason these two gave lesser indications of 
the presence of anomalies was because the focussing effect produces a pattern that covers only a smaller part of the celestial 
sphere, which tends to get lost when averaged over the whole sky. This method is therefore useful to characterize coherent 
signals extended over a large region, such as a spiral pattern, but not if they are concentrated. 

Phase analysis is a relatively new technique, and has consequently not been used to quantify many alternative situations 
to the concordance model. The phases of the spherical harmonic coefficients provide a generic way of looking at correlations in 
harmonic space that could arise from non-stationarity or non-Gaussianity. While this is a potential strength of the approach 
- while phase correlations will not just be useful for identifying anisotropics specific to the Bianchi models, but in theory 
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any isotropy introduced to the CMB - it could also prove a weakness, in that more general methods may lack the power to 
discriminate very specific models. 

The phase correlations identified in our Bianchi maps using this technique were much stronger than we at first expected; 
given the generic nature of the metric it was not expected to yield good results. In addition to this, the strong correlations 
were found to be robust to both rotation and moderate noise. Significant correlations in both twisted and focusing features 
were also identified. However using the same methods on the WMAP 5 year data shows little evidence of non-Gaussianity. 
Given that the diagnostics are identified in harmonic space, it is difficult to say whether any of the anomalies identified this 
way are down to isotropy or homogeneity. 

The analysis of multipole vectors is also a relatively new technique. It has been used to identify non-Gaussianities in the 
WMAP data, and has been particularly successful in identifying anisotropics (i.e. asymmetries and/or preferred directions). 
The multipole vectors are calculated from spherical harmonic coefficients which, as we have already shown, themselves provide 
a very effective way of identifying correlations in Bianchi (and presumably other anisotropic) patterns. The multipole vectors 
must include at least some of the information needed to describe these mode correlations. The advantage of multipole vectors 
over the spherical harmonics themselves, however, is that they give results in real (i.e. pixel) space which is much more 
informative to the user. The results when applied to the Bianchi maps show very strong correlations between the directions 
of the multipole vectors for low often with them entirely overlapping, and hence showing preferred directions. Since these 
vectors would not be aligned in the case of a stationary stochastic field over the sky, these results demonstrate that they are 
sensitive to departures from the standard cosmological model. 

It remains the case that the standard cosmological model is a good fit to a huge range of observational data. Nevertheless, 
it is important that tools are developed that are sufficiently sensitive to hunt efficiently for possible anomalies in the next 
generation of observations. There are many ways that the CMB temperature pattern could be anomalous other than through 
the presence of Bianchi perturbation modes such as those we have studied here. Just as there are many ways a distribution 
can be non-Gaussian, so are there also many ways a fluctuation field can be non-stationary. Testing for departures from 
the standard model will require not one but a battery of statistical techniques each sensitive to particular aspects of the 
distribution. 

This has been a very preliminary analysis, aimed at establishing whether the diagnostics described in this paper are in 
principle capable of uncovering evidence of underlying anomalies in CMB data. Of course these patterns represent somewhat 
extreme departures from the standard framework so it is no real surprise that they register strongly in the descriptors used. 
However, in all cases our analysis has involved only a relatively small number of quantities, so the fact that we see quantifiable 
effects emerging is very encouraging. 
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APPENDIX A: MULTIPOLE VECTORS TERMINOLOGY 

Any homogeneous polynomial Fr of degree i in x, y and z can be written as 

Fr{x, y^z) = \ - {aix + hiy + ciz) • {a2X + 622/ + C2z) • • • {a^x + h^y + c^z) + Sr • Gr. (Al) 

The function Gr is homogeneous of degree ^ — 2 and Sr — x'^ -\- y'^ -\- z'^ . All polynomials Fr^ Gr and Sr and all variables 
are real. If we consider the values taken by the polynomial on the unit sphere where Sr — x'^ -\- y"^ -\- z^ — 1 then the above 
expression reduces to the product of the linear parts together with the term Gr. This can be extended to a complex polynomial 
through 

F(x, y, z) = AL1L2 -"L^ + S'G (A2) 

where Li — aiX + biy + CiZ and x, y and z are complex numbers, while the coefficients of F (ai, bi,Ci, and A) are real. G is a 
homogeneous polynomial of degree £ — 2 and S is given by x^ -\- y'^ -\- z'^ . We are interested in the value of the polynomial F 
on the 2-sphere which, in the complex space, takes the form S = x^ -\- y'^ -\- z'^ =0. By Bezout's theorem, which holds that 
the number of points on two curves is equal to the product of their degrees, there are 2i points in which the complex curve 
F{x,y, z) = intersects the quadratic curve S{x,y, z) = 0, which is topologically a 2-sphere. Since the complex curve F = 
intersects the complex S = ^ in 21 points, the product curves Li also intersect 5* = in the same 21 points, fi = {xi,yi, Zi). 
Moreover, its complex conjugate /* = {x* ^y* , z*) has to lie in the intersection in which both F and S are zero. We thus 
obtain the 2i points of intersection {/i, /i , • • • , Z^, //}. 
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Each pair j/^, //} determines a unique hne Li — aiX + bit/ + CiZ = with real coefficients such that: 

Re I 7 Re I Re p. 

aiXi + biyi + CiZi - 

a,x'r + b,y'r + c,z'r = (A3) 

where the coefficients are normalized to unit length, i.e. af -\- bf -\- cf = 1. Note that the index i has no sum and 1 ^ i ^ i. 
In order to find the multipole vectors, Vi = (a^, bi, c^), of each i we need to find the pairs of fi which lie on the curve F = 
and 5* = in the complex projective plane i.e. finding the roots (a) which satisfy F = on the 2-sphere (5* = 0). What is 
required is to factorize the homogeneous, harmonic polynomial, F into linear factors i.e. such that F is the product of Li 
only. However, this is not possible analytically since they are not linear equations far from the dipole. Fortunately, the curve 
S = can be parametrized as a single variable and the polynomial F as: 

y, z) = F{i(a^ - 1), -2ia, + (A4) 

From Equation IA4I the roots a which satisfy F = and 5* = 0, or the product of Li=0, can be found. Once the roots a 
are found, the pair of {fi, /*} can be expressed in terms of x, y and z. The next step is to find the multipole vectors from 
Equation [A3] by using fi — {x^^,y^^, z^^) + i{x^^ ,y^^ , z^^) or its conjugate /* since the two points give same result. 

We now apply this terminology to the relevant cosmological application, that of temperature patterns on the CMB sky, 
as follows. The multipoles, T^, can be represented by a polynomial F on the 2-sphere (5* = 0), or the product of Li. By 
the given relations x — i{a^ — 1), y — —2ia, and z = + 1, the spherical harmonics can be expanded as a terms. Thus 
the dipole (Ti), quadrupole (T2) and octopole (T3) can be described in terms of a with the aim as coefficients. The aim were 
calculated using Healpix from the maps, and we found the roots a which satisfy equations T^=0 for each i. These roots gave 
a pair of {fi, fl} , therefore from Equation IA3I we obtained solution sets Vi = (a^, bi, Ci) which are multipole vectors for each 
i. 

We now explain this procedure in detail for each of the multipoles considered: 
AO.l Dipole: £=1. 

Tl = ^airnl^lm (A5) 

— (an x-aii y+ —z) (AG) 

For the dipole vector it is not necessary to use the a notation since it is a linear equation in x, y and z. In our polynomial 
notation, the dipole is also represented as 

F = Li = Ai(aix + biy + az). (A7) 
From the normalization, we obtain the multipole vector. 



^1 (\l , -\ ^ciiT, -^==aio), (A8) 



3Ci ' V 3Ci ' V3CT 
and the coefficient. 



where the Ci = (2|aii|^ + a?o)/3 is the angular power spectrum d of the monopole. 



AO. 2 Quadrupole: £^2. 

In this case we need to expand T2 using the a notation since for multipoles from the quadrupole to higher order the equations 
are no longer linear. The quadrupole on the sphere has two multipole vectors, vi and V2, from F = L1L2 = Ai(aix + biy + 
ciz){a2X + b2y + C2z). In order to find them, we transfer the quadrupole expression from spherical harmonics to a notation 
for efficient computing, 

T2 = ^a2ml2m (AlO) 
Re I o- Re\ 4 a/ Im ■ Im\ 3 , / /T^ , r' Re\ 2 

2«20 - ^22 + 2^a2l - 4(a22 - ia2i )a + (v6a2o + 6032 

I A / Im I • Im\ I / 3 Pie n ■ Re / a -i -i \ 

+4(a22 + ^a2l + y -^20 - ^22 - ^10.21 (All) 
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AO. 3 Octopole: i = 3. 

We use the same method as we have done for the quadrupole. The octopole has three multipole vectors, vi, V2 and vs, from 
F = L1L2L3 — Xi(aix + biy + ciz){a2X + b2y + C2z){a3X + bsy + csz), which gives a 6*^ order of equation in a: 

T3 = ^asmlsm (A12) 

m 

= Aea^ + Asa^ + A4a^ + Asa^ + Asa^ + Aia + Ao, (A13) 
in which the coefficients are, 

Ae = 5a3o - VsOa^i + {-bVSai^' + ^afs")^ 
As = -4V30a^r + (-lOVSa^r + QV^aiT)i] 
A4 = 5[3a3o + VsOa^i - {VSa^i + 3V5a^i)i 
As - -20(V3a^r + V^aiT)i 

A2 = = -a;,Ao = 

Each multipole has 2£ roots of a which gives the (fi^ft) pairs, however only i components are used to find the multipoles 
since their conjugators give the same results, as we mentioned earlier. 
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